Dependence of electronic and optical properties on a high-frequency field for carbon nanotubes 
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We study theoretically the electronic structure, transport and optical properties for a zigzag single-wall carbon 
nanotube connected to two normal conductor leads under the irradiation of an external electromagnetic field at 
low temperatures, with particular emphasis on the features of high-frequency response. Using the standard 
nonequilibrium Green's function techniques, we examine the time-averaged density of states, the conductivity, 
the dielectric function and the electron energy loss spectra for the system with photon polarization parallel with 
the tunneling current direction, respectively. Through some numerical examples, it is shown that the density of 
states is strongly dependent on the incident electron energy, the strength and frequency of the applied field. For 
higher electron energies in comparison with lead-nanotube coupling energy, the system conductance decreases 
with increasing the field strength and increases with increasing the field frequency respectively, and shows some 
oscillation structures. Moreover, the optical functions for the system have also a rich structure with the variation 
of field frequency. It may demonstrate that this transport dependence on the external field parameters can be 
used to give the energy spectra information of carbon nanotubes and to detect the high-frequency microwave 
irradiation. 

PACS numbers: 73.23.-b, 73.63.Fg 



I. INTRODUCTION 



Nanoelectronics, or molecular electronics, have been 
proposed as the alternative to silicon in future technical 
applications' and have attracted much interest recently. Es- 
pecially, carbon based nanostructures, such as fullerene, 
graphene and nanotubes, are the most interesting structures 
because of their excellent physical and chemical proper- 
ties. Carbon nanotubes (CNTs) can exhibit metallic or 
semiconducting behavior with different diameter and chi- 
rality, and therefore they will be promising candidates for 
the new carbon-nanotube-based electronic devices, such as 
nanosensors,^ electric batteries,^ field-effect transistors,"''^ 
Coulomb blockade devices, ^'^ and field-emission displays.** 

The quantum transport properties of CNT systems have 
been investigated experimentally by many authors.^ " The 
coherent quantum wire behaviors of an individual single- 
walled carbon nanotube (SWCNT)'^ between two contacts and 
the conductance quantization of multi-walled CNTs'° have 
been observed, respectively. Moreover, Onac et al}^ have 
used an on-chip to detect the high-frequency noise signals 
generated by quantum dot formed in a SWCNT with photon- 
assisted tunneling. This result so far indicates no intrinsic fre- 
quency limitation due to the CNTs themselves. Using the 
long-wavelength approximation, Mele et al}^ have studied 
the coherent one-photon and two-photon electronic excita- 
tions for graphene sheets and SWCNTs. The optical dielec- 
tric function of a finite SWCNT has also been studied early by 
gradient approximation'^ and by first-principle calculations,'^ 
respectively. Furthermore, the conductances for microwave 
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field irradiated multiwall CNT'^ and SWCNT"' connected to 
two leads have been investigated by photon-assisted transport 
model, in which the conductance oscillations evolve toward 
a well-defined step structure and its sensitive dependance on 
the field intensity have been expected.'^ This study has been 
extended'^ to optical selection rule for SWCNT very recently. 
The question thus arises whether SWCNTs could be useful 
as terahertz (THz) detection, which is more important in the 
use of quantum control to nanostructures. To the best of our 
acknowledge, this important case for a driven SWCNT-based 
system has not been reported previously. 

In this paper, we present a theoretical investigation of the 
electronic structure, transport and optical properties for a 
metallic zigzag chiral SWCNT under a THz electromagnetic 
field (EMF) irradiation since metallic CNTs are known to be 
one of the ideal system exhibiting quantum mechanical nature 
of the electrical transport. ' '' The dependence of time-averaged 
electron density of states (DOS), conductance, dielectric func- 
tion (e) and electron energy loss spectrum (EELS) on the ir- 
radiation field strength and frequency are demonstrated for 
the system, respectively. Some different characteristics is ob- 
tained and the results are compared with those for the similar 
physical systems in the previous works. ^•'■'"'^ 

The rest of the paper is organized as follows. The analytical 
expressions of the time-averaged DOS, the conductance and 
the optical functions are calculated starting from the system 
Hamiltonian by NGF approach in Sec. II. Some numerical 
examples and discussions for the results are demonstrated in 
Sec. III. Finally, Sec. IV concludes the paper. 



II. MODEL AND FORMULISM 

The NGF approach has been employed in last decades to 
study transport problems involving variety of interactions be- 
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yond the linear response regime. Meir et al. have derived 
a formula for the current through a region of interacting elec- 
trons using the nonequilibrium Keldysh formulism. Here we 
assume that the external high-frequency field only causes a 
rigid shift in the electron energy spectrum under the widely 
adopted adiabatic approximation. Therefore, the Hamiltonian 
for an irradiated SWCNT embedded between two normal con- 
ductor leads reads 

q,a€L/R k 

+ J] ^'^l-^'kdl^'^k + Vl^j^C^dq^a) 

q,aeL/R,k 

+ ^ ^0 cosia)t)c^Ck', (1) 

k,k' 

where the operator dg „ (dg^a) creates (annihilates) an electron 
with momentum q in mode a in either left (L) or right (R) 
lead, and (ck) creates (annihilates) an electron at the state k 
of the CNT with energy spectrum'^-^" ek-+j[3 + 2 cos(2kyb) + 
4 cos(/:v/^) cos(3fcifl)]'^^, which is similar as that of graphene 
though the boundary conditions of SWCNT are different from 
graphene, especially in the circle direction. Where y is the 
hopping integral of CNT in the tight-banding approximattion, 
a siO.TlA is the half length of in-plane lattice constant with 
b= "V^a, and kx and ky are in-plane wavevectors of Dirac elec- 
tron. The coupling between the electrode leads and the CNT 
with strength Vq^a,k is represented by the third term in Hamil- 
tonian (1), and the last term describes the electron-photon in- 
teraction in the CNT, where Vq is the interaction strength un- 



der dipole approximation and co the incident field frequency. 
Moreover, for simplicity, in the above Hamiltonian we have 
neglected the spin and electron-electron interactions although 
much attention has been paid to the pseudospin-related eff'ects 
in graphenes^' and SWCNTs.^^ 

Now we employ the usually defined retarded and corre- 
lated Green's Function (GF)"^-''^'!^ as G[(t2,ti)=+i0(±t2 + 
hXickitilc^iti)]) and G<(f2,fi)=-!X{Q(fi),c+(f2))> to the 
Hamiltonian. When the unperturbed GF g' ie) of the nanotube 
is obtained, one can obtain GFs of Cie) and G^(e) from the 
Dyson equation G''(e)=[(^''(e))"> - 2;''(e)]-i and the Keldysh 
equation G^(e)=G''(e)S^(e)G''(e), respectively. Further- 
more, the self-energies are E'^(e)=(-i/2)[E^(e) + S*(e)] and 
2<(e)=(-i/2)[/z,(e)2^(e) + /R(e)E*(e)], where the linewidth 
function S°'(e)=27r Y.g,k ^q,a,k^q,a,k^i^ - ^q,a) describcs the in- 
fluence of the leads and fa(e) is the Fermi distribution func- 
tion. In the wide -bandwidth approximation,^^ '*'^^ '^ 2"(e) is 
taken to be independent of the energy and energy levels. 

Furthermore, the GF of a honeycomb lattice (graphene) can 
be written as 



^ ^2 _ 3 _ 2 co^{2kyb) - 4 co%(kyb) co^iZk^a) ' 

For a (n,0) chiral zigzag SWCNT, the periodic boundary 
conditions'^'^^ requires that the transverse wave vector /cj, sup- 
plies kyb=nj/n with good quantum numbers for the subbands 
j=l~2n, and the longitudinal wave vector k^ is confined in 
the first Brillouin zone -n/ {2>a)<kx<nl (3a). Applying Fourier 
transform to g^(e) one obtains the unperturbed GF in the real 
space 



2n ^;r/(3a) 



g\e, I, m) : 



J-./(3a 



dk. 



-3ikja -injm/n 



€^ -3 - 2 cos(2;r j/n) - 4 cos(;r j/n) cos(3kxd) ' 
I 



(3) 



where / and m are respectively the length and width index in 
real space, once the integration over kx is done in Eq. (3), 
one can obtain the analytical expression of Cie) from Dyson 
equation. Therefore, the real part of the zero-bias limit linear 
conductance^^"'^ can be derived as 



^t^G'-(e)S«G''(e), 



(4) 



where J^ix) is the /ith-order Bessel function of the first kind. 
The imaginary part of conductance, 0-2(1^0;^^)^ can be calcu- 
lated from CTiiVo, cj) by Kramers-Kronig transformation'"*''^ 



criiVo,oj) = 1 -I- 



ip r 

^ Jo 



do) 



ll> o-2(Vo, oj ) 
~ 2 ' 



(5) 



complex dielectric function e(Vo,a>)=l+4nio-(VQ,uj)/a) with 
cr(Vo,w)= o-i(Vo,io)+o-2(Vo,co), in principle one can obtain 
all other linear optical properties such as reflectivity, absorp- 
tion spectrum and EELS by -3m[e(Vo, a>)]"' at the long- 
wavelength limit,'-' ''* where 3m represents the imaginary part 
of the quantity. The effective electron DOS will be'^ 



Do(£ - fiHcj), 



(6) 



In the tight-banding approximation, the DOS of a (n, 0) CNT 
in the absence of perturbing potential reads^^ 

Z)o(e) = — 3m^jXe-i-«0"')x|8[l-i-cos(2;r77n)] 



7=1 



[(£ + iO^f -2cos(2nj/n) - 3f] 



(7) 



where P denotes the principal value of the integral. From the 



It should be pointed out that Eq. (6) can be interpreted as 
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following. Photon absorption (fi>0) or emission (fi<0) can be 
viewed as creating an effective electron DOS at sideband en- 
ergies e^-jJiu with a probability \Jfi{Volfio})^. One notes that 
the multi-photon processes will be suppressed more or less 
in the cases of Vo/(^w)»l, and the enough strong EM field 
(Vo/(^^<j)»l cases) will split sharp singularities in the elec- 
tron spectrum except for its aff'ection on the electron distri- 
bution of the CNT.-^^ Therefore, in this work we concentrate 
on the cases of the relatively high-frequency response under 
moderate irradiation field strength for the system. 



III. RESULTS AND DISCUSSIONS 

In the following, we present some numerical examples of 
the calculated normaUzed DOS, linear conductance and opti- 
cal properties for an external field irradiated (12,0) SWCNT. 
The tight-banding approximation hopping integral 7=2.75 eV 
is taken to be the unit of energy, the temperature^"* and the ef- 
fective 'speed of light' are A:sT=0.01 and vf=8.88x10^ m/s, 
respectively. In our numerical calculation, the geometry in- 
dex / and m in Eq. (3) are selected as 12 and 24 for a 
SWCNT system, that is, the diameter and length of the CNT 
is about 1 nm and 2.6 nm, respectively. In the wide-band 
approximation,'^''*'^**'''' we assume that 2^=2^=2=0.0017 as 
the CNT-lead coupling parameter for simplicity. The above 
selected parameters for the system are reachable in the present 
experiments.^"" 

Figure 1 shows the calculated time-averaged DOS (in arbi- 
trary units) and linear conductance ai (in units of <To=2e^//j) 
as a function of the incident electron energy e (in units of 
y) with a fixed field frequency v=(y/(27r)=33.25 THz for dif- 
ferent field strength Vq. In the absence of irradiation field, 
as shown by the soUd Une in Fig. 1(a), the DOS presents a 
series of original peaks at energies ±0.40, ±0.50, ±0.75 and 
±0.95, respectively. A wide van Hove pseudogap is observed 
wherever the DOS is not zero. However, when an irradia- 
tion field of strength yo=l.l eV with photon energy hv=0.Q5y 
is applied, as shown by the dashed Une in Fig. 1(a), ex- 
cept for a series of resonance peaks in low energy range of 
-0.407<e<0.407, there are some resonance peaks at ±0.55, 
±0.60, ±0.65, ±0.70, ±0.80, ±0.85 and ±0.957, respectively. 
These additional peaks are mainly attributed to the processes 
of 1-photon and 3-photon processes. Moreover, in the case of 
field strength yo=2.20 eV [see the doted Une in Fig. 1(a)] with 
the same field frequency, due to the 2-photon, 4-photon and 
6-photon processes, the system DOS presents a much richer 
structure except for a remarkably changed van Hove pseudo- 
gap. This result is quantitatively diflrerent from that in Ref. 16 
(see, e.g.. Fig. 2 in Ref. 16) where no quantitative explana- 
tion about the electronic structures was given. The sensitive 
dependence of DOS on the irradiation field strength can be 
physically interpreted as following. With a certain photon fre- 
quency, the variation of the field strength should be followed 
by diff'erent photon sideband processes as can be verified from 
Eq. (6). Furthermore, the electronic structures of an irradi- 
ated (10,0) zigzag SWCNT (semiconducting type) have also 
been plotted with the same parameters (not shown here), the 




FIG. 1: The dependence of (a) time-averaged DOS and (b) con- 
ductance on the incident electron energy e with a field frequency 
v=33.25 THz for three different field strengths, where solid line 
[down-shifted 0.1 for comparison in (a)] for Vo=0, dashed line for 
Vo=l-l eV and dotted line [up-shifted 0.1 for comparison in (b)] for 
Vo=2.2 eV, respectively. 



results are similar except for a wider van Hove pseudogap is 
observed. 

The corresponding real part conductance (r\ to Fig. 1(a) 
is illustrated in Fig. 1(b). In the absence of external filed, 
as shown by the solid line in Fig. 1(b), cri presents a peak 
at the four original resonance states and cri is nearly in the 
pseudogap as expected. Wherever, under the irradiation of a 
field with strength of yo=1.10 eV, as shown by the dashed 
Une in Fig. 1(b), the system conductance presents a series of 
additional peaks at ±0.55, ±0.60, ±0.65, ±0.70, ±0.75 and 
±0.857, respectively. When the field strength increases to 
2.20 eV, one can find several additional conductance peaks at 
±0.15, ±0.20, ±0.25 and ±0.307 i^i the low energy range [as 
shown by the doted line in Fig. 1(b)] owing to multi -photon 
processes, while the conductance peaks in the higher energy 
range are unchanged. 

In what follows, we show the influences of the external field 
parameters on the electronic structure and transport of the sys- 
tem for a fixed incident electron energy 6=87/1 1 (2.0 eV). Fig. 
2 depicts the time-averaged DOS and <ti as a function of the 
irradiation field strength (~ Vq). In the case of a relative low 
field frequency v=33.25 THz, as shown by the solid line in 
Fig. 2(a), the DOS decreases slowly with a oscillating period- 
icity about 0.5 eV other than a step at 0.50<yo<l-0 eV. When 
the field frequency increases to v-49.9 THz (dashed line) and 
v=66.5 THz (dotted line) as shown in Fig. 2(a), the oscillating 
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FIG. 2: The dependence of (a) time-averaged DOS and (b) conduc- 
tance on the field strength with a fixed energy e=2.0 eV for three 
different field frequencies, where solid line for v=33.25 THz, dashed 
Une for v=49.9 THz and dotted Une for v=66.5 THz, respectively. 



FIG. 3: The dependence of (a) time-averaged DOS and (b) con- 
ductance on the field frequency with a fixed incident electrons en- 
ergy e=2.0 eV for three different field strengths, where solid line for 
Vo=0.55 eV, dashed line for Vo=l l eV and dotted line [up-shifted 
0.2 for comparison in (a)] for Vo=2.2 eV, respectively. 



periodicity seems to become larger except for a peak at 0.9 eV 
and 1 .0 eV, respectively. 

Figure 2(b) illustrates the behavior of cri as a function of 
the field strength for three different field frequencies. In the 
case of low frequency v=33.25 THz, as shown in the solid 
line, (Ti presents a downgrade step-like structure with oscil- 
lating periodicity about 0.5 eV similar to that of the DOS [see 
Fig. 2(a)]. When the field frequency increases to v=49.9 
THz or v=66.5 THz, as shown in the dashed line or doted 
line of Fig. 2(b), the oscillating periodicity becomes larger. 
As mentioned above, due to the case of higher frequency for 
the system we have considered, the results here are quanti- 
tatively different from those in Ref. 16 (see, e.g.. Figs. 4 
and 6 in Ref. 16) where well-defined up-going step structures 
for the conductance have been predicted only for very low 
(y<3.3 THz) field frequencies with zero incident electrons en- 
ergy and temperature. It should be pointed out that the results 
here are qualitatively different those in Ref. 15 (see, e.g., Figs. 
2 and 3 in Ref. 15) where a nearly universal down-going nor- 
malized resistance with the microwave power and frequency- 
irrelevant normalized resistance at liquid helium temperature 
have been observed. Moreover, the electronic structures and 
conductance of a longer zigzag SWCNT have been plotted 
(not shown here) with the same all other parameters, the re- 
sults are qualitatively similar as expected. 

Figure 3 demonstrates the calculated time-averaged DOS 
and conductance cr\ as a function of the irradiation field fre- 
quency V (in units of THz) with a fixed incident electron en- 



ergy 6=2.0 eV (same as in Figure 2) for three different field 
strengths. As shown in Fig. 3(a), one notices that in all case 
of field strength there is a series of peaks at frequencies v=10, 
12, 14, 16.5, 20, 22.5, 26, 30, 36, 45.5 and 60 THz, respec- 
tively. Owing to the difference of the field strengths, these 
peaks should result from different photons processes as men- 
tioned above and the characteristics will be manifested in the 
conductance more or less. The behavior of the time-averaged 
normalized conductance iti versus the field frequency for dif- 
ferent strengths is illustrated in Fig. 3(b). In the case of 
smaller strength Vo=0.55 eV (see the solid line), cri presents 
a rapidly oscillatory up-going structure with the increase of 
field frequency, and reaching almost 2.0cro as field frequency 
reaches 36.5 THz. However, when Vq increases to 1.10 eV, as 
shown the dashed line in Fig. 3(b), cr\ up-goes slowly, only 
two obvious upgrade step-like structures at 20 and 30 THz, re- 
spectively. Furthermore, even a much more slowly oscillatory 
structure is shown by the doted line in Fig. 3(b) for the case 
of lager Vo=2.20 eV. It should point out that the series of res- 
onance peaks are due to multi-photon processes. The above 
results are quantitatively different from those in Ref. 16 (see, 
e.g.. Fig. 5 in Ref. 16) where a nearly independence of field 
frequency in the range of 0<v<0.0057 and a higher frequency 
0.0057<y<0.05y oscillatory response of the conductance have 
been predicted. 

Next, we present numerically the dielectric function s(v) 
and EELS versus the irradiation field frequency with a fixed 
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FIG. 4: The dependence of (a) dielectric function e(v) and (b) elec- 
tron energy loss spectra (EELS) on the field frequency with a fixed 
energy 6=2.75 meV for two different field strengths, where solid line 
[up-shifted 0.5 for comparison in (a) and (b)] for Vo=l.l eV and 
dashed line for Vo=2.2 eV, respectively. 

incident electron energy (selected as 2.75 meV for compar- 
ison with the results in Refs. 14-17) for two different field 
strengths, and further discuss the basic optical properties '^'"'^ 
of the system. As shown in Fig. 4(a), e(v) in dashed line, one 
notices a broad optical absorption frequency strip at the height 
of 1 .0 and 2.0 in the ranges of 16.5<v<23.5 THz and 24<v<30 
THz respectively, with down-going oscillation structures in 
the range of 30<v<50 THz for the case of yo=2.2 eV. Wher- 
ever, a peak of near 2.0 height in the vicinity of 13.5 THz 
with down-going oscillation structures in the higher frequency 
range for yo=l.l eV (see the solid line). The characteristics 
of the s{v) are different from the results in Refs. 13, 14 and 
17 since different parameters we have chosen, although the 
principal features of Fig. 8 in Ref. 14 have been revealed. 
Corresponding to the s(v) as show in Fig. 4(a), the EELS 
is illustrated in Fig. 4(b). For the case of yo=2.2 eV EELS 
presents a richer structure as expected since it is related to e(v) 
by -3m[e(v)]"' . Except for the two broad peaks in the neigh- 
bor of 10 and 12 THz for y()=l.l and 2.2 eV, one notes a op- 
tical absorption strip at the height of 1.5 and 0.8 in the ranges 
of 13.5<v<16 THz and 17<v<23 THz respectively, with os- 



cillation structures around 0.5 in the range of 24<v<50 THz 
for the case of yo=2.2 eV [see the dashed line in Fig. 4(b)]. 
Since different nanotubes always have a dominating plasmon 
mode'^ at ~6 eV, therefore, all the above features should rely 
on multi-photon processes. In short, all the electron DOS, 
conductance, dielectric function and EELS of CNTs system 
are sensitive to the strength and the frequency of the external 
field, which may be applied to the detection of the external 
THz field and the energy spectra information of carbon nan- 
otubes. 



IV. CONCLUSION 

In summary, using the NGF method, we have investigated 
theoretically the electronic structure, transport and optical 
properties of a metallic-typed zigzag single-wall carbon nan- 
otube connected with two normal leads under the irradiation 
of an external high-frequency (THz) electromagnetic field at 
low temperatures. Using the standard nonequilibrium Green's 
function techniques, we examine the time-averaged density 
of states, conductance, dielectric function and electron energy 
loss spectrum for the system with photon polarization parallel 
with the tunneling current direction, respectively. It is demon- 
strated that, by analyzing some numerical examples, the den- 
sity of states shows a strong dependance on the incident elec- 
tron energy and external field parameters. For high irradiation 
energies in comparison with CNT-lead coupling, the system 
conductance decreases with the field strength while increases 
with the field frequency followed by some oscillation struc- 
tures. For the lower incident electron energy, the dielectric 
function and electron energy loss spectrum are consistent with 
the results in Refs. 13, 14 and 17 due to multi-photon pro- 
cesses. In all cases the system transport and optical properties 
are sensitive to the parameters of the external field, which may 
be utilized to the detection of high-frequency irradiation and 
the energy spectra information of carbon nanotubes. However, 
the experimental observation for these effects and further the- 
oretical investigations on the system with impurity, spin or 
electron-electron interactions and the armchair chiral carbon 
nanotubes are worthy to be carried out. 
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